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Abstract 

Deconfinement and screening of higher-representation sources in finite-temperature 
SU (2) lattice gauge theory is investigated by both analytical and numerical means. 
The effective Polyakov-line action at strong coupling is simulated by an efficient 
cluster-updating Monte Carlo algorithm for the case of d = 4 dimensions. The re- 
sults compare very favourably with an improved mean-field solution. The limit d^oo 
of the SU(2) theory is shown to be highly singular as far as critical behaviour is 
concerned. In that limit the leading amplitudes of higher representation Polyakov 
lines vanish at strong coupling, and subleading exponents become dominant. Each of 
the higher-representation sources then effectively carry with them their own critical 
exponents. 
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For non-Abelian gauge theories one can define a static potential between matter sources 
transforming as arbitrary irreducible representations of the gauge group, from now on taken 
to be SU(N). This static potential between two infinitely heavy sources is believed to 
depend crucially on the manner in which the chosen representation behaves under trans- 
formations restricted to the center Z(N) of the gauge group. Representations that are 
insensitive to Z(N) transformations should yield a screened potential, while those sensitive 
to these transformations should yield a confining potential. This is the standard picture of 
confinement and screening in non-Abelian gauge theories, dating back almost twenty years 
(see, e.g., ref. 0). 

With numerical simulations these ideas can be tested in the context of lattice gauge the- 
ories. Indeed, one can go further and investigate in more details the dynamics behind these 
screening and confining mechanisms. In particular, one can measure the relevant distance 
scales, those that separate "short distances" (essentially perturbative physics, on account 
of asymptotic freedom) from "long distances" (a confining string between two irreducible 
sources transforming non-trivially under Z(N), a screened potential between those trans- 
forming trivially). The dynamics of the intermediate region has been found to be, in many 
respects, quite rich. It has been observed that even representations transforming trivially 
under the center group could feel a linearly rising potential (with a slope different from the 
string tension of the fundamental representation) at intermediate distances @ . At a certain 
range of distance scales all representations appear to carry with them their own dynam- 
ics. In the limit of an infinite number of colours, factorization is sufficient to show that all 
irreducible representations are confined by a linearly rising potential if the fundamental rep- 
resentation is, with string tensions that depend on the representations ||. Essentially, the 
intermediate distance region in which a non-zero string tension exists for all representations 
grows with N, the number of colours, reaching infinity as iV — > oo. 

These results ought to have some bearing on the physics of finite-temperature gauge 
theories as well. Indeed, some very simple numerical simulations have shown that just as, 
for example, adjoint sources may feel a linearly rising potential at intermediate distances, 
such sources appear to "deconfine" at precisely the same critical temperature T c at which 
fundamental sources deconfine |4[].[] Since, for finite N, adjoint sources are not genuinely 
confined out to arbitrarily large distances this is of course only to be understood at the 
qualitative level. 

For gauge theories with continuous deconfinement phase transitions one obvious question 
concerns the critical behavior of the Polyakov lines corresponding to sources of arbitrary rep- 
resentations. Again, representations transforming trivially under Z(N) should not be able 
to serve as order parameters for the phase transitions, while those sensitive to Z(N) should, 
since it is precisely a global Z(N) symmetry which is broken at the deconfinement phase 
transition JJ. Universality arguments would, for such continuous phase transitions, place 
the finite-temperature SU(N) gauge theory in the universality class of globally Z(N) in- 
variant spin systems with short-range interactions JFJ. The analogue of the spin operator 
would be the Polyakov line in the fundamental representation. What about Polyakov lines 
corresponding to traces taken in higher representations of SU(N)? There appears to be 
no room for independent exponents for these higher representations from the spin-system 

■"■Even representations transforming non-trivially under Z(N) could in principle deconfine at temperatures 
different from the temperature at which the fundamental representation deconfines. This has not been 
observed, however H. 
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point of view, since there is simply no obvious analogue of the "internal" SU(N) degrees of 
freedom in the Z(N) fixed-point language. This would indicate that all higher representa- 
tions sensitive to the center group should be equivalent order parameters with exactly the 
same critical behavior as the fundamental representation. Physically, this would also be in 
agreement with the screening argument according to which all higher representations sensi- 
tive to Z(N) eventually, for large enough distances, should be colour screened as far down 
as possible (i.e. down to the charge in the fundamental representation). On the other hand, 
since strings with representation dependent string tensions do form at some intermediate 
distance scales for all representations, it is not quite obvious how much of the screening 
mechanism will be observable at standard present-day lattice sizes used for Monte Carlo 
simulations. The clear change in behavior of the adjoint Polyakov line at the deconfinement 
"phase transition" (there is of course no genuine transition in a finite volume) as measured 
on small lattices j|] is already one indication that there may be difficulties with numerical 
investigations of this problem. With the same level of statistics (and the same lattice sizes) 
that are used routinely to confirm the universality arguments based on the fundamental 
Polyakov line, a surprisingly different behaviour was found for the higher representations of 
SU{2) lattice gauge theory in (3+1) dimensions || ||. (For the continuous deconfinement 
transitions of SU (2) and SU (3) lattice gauge theories in (2+1) dimensions, see ref. |TD|, [TT|]). 



These numerical simulations indicated that sources of higher representations that were sen- 
sitive to Z(N) would correspond to different magnetization exponents, one exponent for 



each representation. But these results could all be criticized || [L0| on the grounds that 
they also seemed to indicate critical behavior for Polyakov lines that simply could not be 
order parameters for the transition, those of transforming trivially under Z(N). Indeed, in 
a Monte Carlo study of (3+l)-dimensional SU{2) lattice gauge by Kiskis fl2| the expected 
behavior (adjoint Polyakov line non- vanishing across the transition, the isospin 3/2 Polyakov 
line behaving like the fundamental) was eventually extracted very close to the finite-volume 
"critical point" T c . Some of the numerical difficulties involved are discussed in ref. [|Hi|j.p| 

It is important to realize that all of these issues can be addressed even in the strong- 
coupling region of the lattice theory. In fact, in this regime the universality arguments are 
even strengthened (since the effective Polyakov-line interactions can be shown explicitly to 
be short-ranged |i~5], |l6j), and the question of the critical behavior of higher- representation 
sources near the phase transition point is as meaningful in the strong-coupling regime as 
near the continuum limit. The advantage of going to the strong-coupling regime is of course 
that the question here can be studied in a much simplified setting which still captures all 
the essentials. The leading-order effective Polyakov-line action reads, with Tr{W indicating 
the trace in the fundamental representation ||T6f : 

S eff [W] = \jY,{ Tr iW{x)Tr l W\x+ ]) + Tr 1 W\x)Tr 1 W( y x + j)\ . (1) 

Here the sum on j runs over nearest neighbours. The effective coupling J is related to 
the gauge coupling g and N T , the number of time-like links in the compactified temporal 
direction. To lowest order, for £77(2), it is 

2 For a recent discussion of the value of the adjoint Polyakov line at the phase transition point, see also 
ref. M. 
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with I n indicating the nth order modified Bessel function. For SU(2) we will use a notation 
in which Tr n W means the trace taken in the representation of isospin n/2. Tr-^W is thus 
the trace in the adjoint representation, etc. Higher orders in the expansion (1) (and cor- 
rections to the effective coupling (2)) can be computed in a systematic expansion fl7j , but 
we will not need these corrections for the present purpose. The effective action (1) becomes 
asymptotically exact in the strong coupling limit. 

For what follows, it is useful to write the effective Polyakov-line action (1) for SU(2) in 
terms of a new variable = \Tr\W(x). The partition function then takes the following 
form: 

A.7£©(§)©(§ + I) + E^© e ] , (3) 

§,] § 

with a local potential V[& 2 } = |ln [1 - <3>(x) 2 ]. 

There are two simple limiting cases in which the effective Polyakov-line action (1) can be 



Z = / [[© exp 



solved exactly. One is the large- N limit , where the deconfinement phase transition turns 
out to be of first order (in agreement with large- N reduction arguments based directly on 
the full gauge theory Jl9|]), and where universality arguments hence cannot be addressed.^ 
The other exactly solvable case is the mean- field limit in which d — > oo, with d being the 
number of spatial dimensions ||. In this limit one finds a genuine second-order critical 
point for the gauge group SU{2) at a critical coupling J c — > as d — > oo. For J > J c all 
higher-representation expectation values (Tr n W) are non-zero ||: 

< Tr ^> = (" +1 >W' (4) 

with a = 2dJ(TriW) being the self-consistent mean-field solution for the fundamental 
representation. Surprisingly, they all display non-trivial critical behavior close to J c : 

(Tr n W) ~ (J-J c f" , (5) 

where f3 n = n/2. For the fundamental representation this just corresponds to the mean- field 
Ising magnetization exponent f3\ — /3 = 1/2, in complete agreement with the universality 
arguments. For the higher representations this new critical behavior is highly unexpected. 
In this unphysical but exactly solvable limit all standard screening arguments appear to 
break down, and we are seeing new behaviour which is not predicted by universality.^ 

This d^oo result (5) is disturbing on several counts. It is normally assumed that the 
relevant Z(2) spin system universality class to which the SU{2) finite-T phase transition 
should belong (if continuous) would display "classical" mean field exponents all the way 
down from d — oo to the upper critical dimension d u (in this case with d u = 4, the critical 
behaviour being modified by logarithmic corrections just at d — d u ). At a first glance 
this might seem to indicate that the non-trivial behaviour (5) for all representations should 

3 Still, the large- TV solution does display a number of interesting features such as the simultaneous de- 
confinement of all higher representations at the transition temperature T c , independently of whether these 
transformations transform trivially under the center symmetry (in this case U(l)) or not. All representa- 
tions of the Polyakov line are hence equally good order parameters in this special case, and all display a 
discontinuous jump at the phase transition. 

4 One also finds separate exponents S n =3/n for the behaviour of (Tr n W) ~ /i 1 / 5 '* at J = J c in a small 
magnetic field coupled to Tr{W . 
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remain valid for all d > 4 (d = 4 just being the limiting case), with non-trivial critical scaling 
even for representations of integer isospin, and with new critical exponents for all isospin 
half-integer representations as well. The first conclusion simply cannot be correct (and we 
will show explicitly below why the argument is invalid), because at strong coupling one can 
compute, for example, the adjoint Polyakov line and see that to first non-trivial order in 
1/g 2 it is non-zero. This result should be valid at least up to the phase transition point, if 
this critical point lies sufficiently deep inside the strong-coupling region [fbl. What about 



the odd-n representations? Could it be that they display new non-trivial critical behaviour 
of the kind (5) even at d—41 To investigate this question we have first performed a Monte 
Carlo simulation of the effective Polyakov-line action for SU(2) in d—A spatial dimensions. 
Again, the advantages of using directly this effective action instead of the full SU(2) lattice 
gauge theory are enormous from a numerical point of view (extracting the critical indices by 
conventional means is notoriously difficult for lattice gauge theories due to the large lattices 
and high statistics required), and at the same time the question we are addressing is as 
urgent in the context of the effective action (1) as in the full SU{2) gauge theory. 

We have simulated the model numerically using a cluster algorithm similar to that 
proposed in ref. [ 20[| . The algorithm consists of two parts. First there is a single cluster 



update |2T] of the sign of the field $. The delete probability is given by 

p d = exp(-4J($(a;)$(a; + j) + \$(x)$(x + (6) 

Then, since the cluster update is not ergodic, we have supplemented it with a standard 
Metropolis update that allows changes in the absolute value of the field $(x). The ratio of 
the number of single cluster updates to the number of Metropolis sweeps is a free parameter 
of the algorithm. We have fixed it by the following rule of thumb: for every Metropolis sweep 
one performs approximately (lattice volume) /(average cluster size) single cluster updates. 
The step size of the Metropolis algorithm was chosen such that the acceptance rate was 
approximately 1/2. Ben-Av et al. implemented such a combination of the single cluster 
update and a local heat-bath update for the N T = 1 finite temperature SU(2) gauge theory 
in 3 + 1 dimensions. Critical slowing down was drastically reduced compared to the local 
update procedure. 

We have determined the critical coupling of the model using the fourth order cumulant 
JJ\ — 1 — (m 4 )/(3(m 2 ) 2 ) PBl , where m is the single lattice average over $, i.e. the Polyakov 
line in the fundamental representation. We have first simulated lattices of sizes L = 4, 6, 8, 12 
and 16 at 4J = 0.55, which was our guess for the critical coupling after some preliminary 
simulations with low statistics. We have then measured 20,000 times, with one measurement 
being performed after two Metropolis sweeps plus the corresponding number of single cluster 
updates. The integrated autocorrelation time of m was r m « 1.7 in units of measurement 
steps for all our lattices sizes. We have subsequently computed the fourth order cumulant 
Ui in a neighbourhood of the actual simulation coupling, by reweighting (m 2 ) and (m 4 ) to 
the correct Boltzmann weight 

(m»>(J) = E * mreXp(( - J + Jo ^ } , (7) 
Eiexp((-J + J )S) 

where S is defined by S = JS + ^2 X V. The resulting curves are plotted in Fig. 1. The 
crossings of the cumulant provide estimates of the critical coupling J c . The error should 
vanish like L x l v . 
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Crossings between lines corresponding to L = 4 and 6, 6 and 8, 8 and 12, and, finally, 
between L = 12 and 16 occur at 4J cross = 0.5507(8), 0.5509(5), 0.5511(3), and 0.5507(2), 
respectively. We take the crossings of the Binder cumulant for L = 12 and L = 16 lattices 
as our best estimate of the critical coupling, i.e, 4J C = 0.5507(2). The limiting value 
of this crossing is fixed within a given universality class, in this case expected to be the 



one of the 4-d Ising model. Brezin and Zinn- Justin p5[ have argued that the effective 
potential for dimensions d > 4 at the critical point is given by just a 4 term. It follows 
that the fourth-order cumulant should take the value U — 1 — r2/3 = 0.27052... with 
T2 = r(l/4) 4 /(87r 2 ) = 2.1884396... at the critical point. The values we find for the fourth- 
order cumulant U x itself at the crossings, 0.367(10), 0.372(15), 0.378(14) and 0.350(20) 
are considerably larger than this theoretical prediction, and we see no significant trend 
towards smaller values with increasing lattice sizes. Kim and Patrascioiu [24] find values 
for the fourth order cumulant of the the 4-d Ising model on similar and larger lattices that 



are consistent with our numerical result. Also Binder et al. |26| estimate a value of the 
critical cumulant of the 5-d Ising model that is consistent with the above numbers. In their 
theoretical discussion Binder et al. allow for a finite mass term in the effective potential at 
the critical point. 

Other universal quantities can be extracted from the fourth-order Binder cumulant. Its 
derivative should scale like dU/dJ ~ L x l v at the critical coupling. A fit of the data according 
to this equation leads to v = 0.490(6) , which is consistent within two standard deviations 
with the 4-d Ising value v = 0.5. 

Next, consider the fourth-order cumulant U n for higher representations. Fig. 2 shows 
the results for n — 2. With increasing lattice sizes the cumulant converges toward 2/3 
for couplings both below and above the critical point. The value £7 = 2/3 signals a finite 
expectation value of the observable. We hence clearly see that the Polyakov line in the 
adjoint representation is not an order parameter. It stays finite in both phases. The fourth- 
order cumulant for the n = 3 representation takes values close to 2/3 in the broken phase 
and values close to zero in the high temperature phase; it behaves as an order parameter. 
But the curves do not display crossings close to the critical coupling predicted from the 
cumulant for the fundamental representation. The curve for L = 16 comes close to that of 
the fundamental representation, so we might expect that for still larger lattices the cumulant 
for the n = 3 representation converges towards the fundamental ones, and that the crossings 
can then be observed at the critical coupling J c . The data for the higher representation 
were so affected by errors that no reliable results for the cumulants could be extracted. 

It is instructive to look at the finite-size behaviour of the Polyakov lines of different 
representations at the critical point. Values for n = 1,2,3,4 and 5 are given in table 1. 
The numbers are obtained from reweighting the 4J = 0.55 simulations. Asymptotically 
the value of the odd representations should scale down with increasing lattice size like 
(Tr n W) ~ L~PI V = L^ 1 , while the values for the even representations should converge to 
a finite value. It appears that the values for the even representations indeed stabilize at a 
(very small) finite value. The behaviour of the magnetization of the odd representations are 
best visualized by looking at (Tr n W)L, which in the asymptotic scaling regime should give 
a constant. For n=l this behaviour is nicely seen in the data, for n = 3 there appears to be 
a stabilization for the largest lattice sizes, while for n = 5 no sign of stabilization is visible 
on the lattices we have considered. It thus appears that the higher the representation, the 
larger lattices are required to see the correct infinite- volume behaviour. Physically, this also 
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makes sense in the light of screening considerations. The higher the representation, the 
more screening is required to reproduce the behaviour of the fundamental representation, 
and the larger distances one needs to probe in order to see this. 

We have also simulated the model for various J > J c on lattices of sizes up to L = 
16. Here we have typically performed 10, 000 measurements per simulation, the aim being 
an approximate determination of the critical exponents (3 n directly from the Monte Carlo 
measurements of the different representations. We have not attempted to fit our data to an 
ansatz including corrections to scaling, but have instead defined a J-dependent "effective" 
exponent /3%f by 

The derivative of (Tr n W) with respect to J was computed from the relation 

^j{Tr n W) = ((Tr n W) ■ S) - (Tr n W)(S) . (9) 

The final results are presented in fig. 4 for the odd representations, and in fig. 5 for the 
even ones. We have carefully checked the dependence of the result on the lattice size and 
included only values which were consistent on the two largest lattice sizes considered. It 
turn out that for the coupling close to J c even the L — 16 lattice was not sufficient to give 
a stable result for the n = 5 representation.^ 

Figs. 4 and 5 demonstrate fairly convincingly that the odd representations converge 
toward an effective = 1/2 independent of n, while the even representations converge 
toward j3^/^ = (as expected if these representations remain finite at J c ). But the plots also 
reveal an interesting phenomenon for larger values of (J — J c )/J c : the effective J-dependent 
exponents (5^^ quickly reach a regime of couplings where they are essentially equally spaced, 
growing linearly with n. Although they never actually reach the mean-field prediction (5), 
they get quite close, and they certainly obey the rule (3^^ ~ n ■ (3^ to surprisingly high 
accuracy. This is just as for the original observations in the full (3 + l)-dimensional SU(2) 
gauge theory || ||. It appears that this approximate linear relation between the /3 n 's, when 
measured not too close to the critical point, can be viewed as the "remnant" of the d = oo 
solution. It is then only very close to the critical point the behaviour changes, and the single 
critical exponent (3 emerges for the odd representations, while the even representations run 
smoothly across the transition point. We can estimate this narrow window in the original 
gauge coupling 4/ g 2 by using the relation (2). In the case of N T = 2 the transition occurs 
at 4/g 2 = 1.6424(4). In order to obtain j3 e f f < 0.625 (i.e. 25% above the correct value 
j3 = 0.5) for the n = 3 representation we would have to take 4/g 2 < 1.66. 

While these results may have clarified the situation in the d = 4 theory, we are still left 
with the surprising d = oo results where mean field theory is believed to be exact. How can 
they be explained? Consider the representation of the effective Polyakov-line action given 
in eq. (3). This is a Z(2)-invariant effective scalar field theory in d dimensions, as expected 
on general grounds. But it is a very particular effective scalar theory, one that embodies 
the underlying SU(2) structure (in the restrictions on the integration interval of and 
in the very special form of the local potential V[$ 2 ], which reflects the Haar measure for 

sum 

5 The curves plotted in these two figures are improved mean-field predictions. See below. 
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Since we at this point wish to focus on the d = oo results, we can restrict ourselves to 



"classical" mean-field considerations. It is instructive |1(J to generalize the partition function 
above to an arbitrary local potential V^[$ 2 ] and relax the limitation on the integrations over 
to be in the interval [—1,1]. The d = oo solution is then found by considering the 
single-site partition function 



/oo 
[[©]exp [C© + V[© e ]] 
-oo 



(10) 



where v = 4<iJ($) will be determined by the self-consistency solution. Clearly, for n being 
any non-negative integer, ($ 2n+1 ) = unless the global Z(2) symmetry is spontaneously 
broken. Call the critical coupling at which this occurs J c . If the phase transition is con- 
tinuous, (<&) will be small just above J c , and it is meaningful to expand in v (no matter 
how large d is taken, once fixed). The result is, for the expectation values of the first two 
non-trivial mean-field moments of $ |1T0 



<$ 2 > 

($3) 



<$ 2 >o+2 
(<*>V + 7 



($ 4 >o 



<$ 2 >< 



v z + ... 



;($ 6 ) -(<l> 2 )o($ 4 )c 



V 3 + 



where the subscript "0" indicates the (constant) expectation value in the unbroken phase 
J < J c . Higher moments can be worked out analogously, by expanding both the partition 
function Z$s an d the unweighted averages in powers of v. Using the recursion relation 
Xn+i = XnXi ~ Xn-i f° r SU (2) characters, we find the general d = oo predictions JLO 



(Tr 2 W) 
(Tr 3 W) 



4($ 2 ) - 1 
A 2 + B 2 v 2 + ... 

8($ 4 )o-4($ 2 ) ] v + A 
A 3 v + B 3 v 3 + ... , 



($ 4 )o-(($ 2 )o)' 



v 2 + 



^$ 6 )o-($ 2 )o($ 4 )o + i($ 2 )o' 



v 3 + 



:i2) 



where A 2 ,B 2 ,A 3 and B 3 are (non-universal) constants. This shows the behaviour expected 
from universality arguments. The adjoint Polyakov line will remain non-vanishing across 
the phase transition at J c (and is hence not an order parameter), and the isospin-3/2 rep- 
resentation scales near J c as v, i.e., as the fundamental representation. But if we take the 
particular potential \^[$ 2 ] of eq. (3), and restrict the integration over $ to the interval 
[—1,1], then devious cancellations occur. One finds ($ 2 )o = 1/4 and ($ 4 )o = 1/8, leading 
to 



(Tr 2 W) 
(Tr 3 W) 



2d 2 J 2 ($) 2 + ... 



v 3 + 



(13) 



It is thus suddenly the non-leading terms in the general expansion of the Polyakov lines that 
become important, due to the amplitudes of the leading terms vanishing in this limit. The 
cancellations required for this phenomenon are actually simple consequences of the orthog- 
onality relations for SU(2) characters, as follows if one performs the mean field calculation 
directly in SU(2) language ||. They occur similarly for all higher representations, leading, 
of course, eventually to the general d = oo solution (5). 
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It is interesting to compare this result with a general renormalization-group analysis by 
Kiskis |T^]. Not being restricted to d — oo, the results of Kiskis can be summarized by 

(Tr 2 W) = a 2 + b 2 t 1 ~ a + ... 

(Tr 3 W) = a 3 t? + b 3 t 1+ P + ... , (14) 

where a is the usual specific heat critical exponent, and t is the reduced temperature near 
the critical point at T c . In the d = oo limit, these are precisely of the form (12) above 
provided we make the identifications 

I- a = 2(3, l+/3 = 3/3 (15) 

(assuming that the coefficients B 2 ,B 3 ,b 2 and 63 are all non-zero). Solving these equations, 
we find (3 = 1/2 and a = 0, the d = 00 Ising model exponents. With the above qualification, 
the identities (15) appear to be new scaling relations for the Z(2) fixed point at d = 00, 
imposed by the combined restrictions of mean field theory and the renormalization group. 

We are now in a better position to understand the d = 00 results. As shown above, 
the appearance of new exponents for each of the odd-n representations in the limit d = oo 
is due to very delicate cancellations that make the amplitudes of the leading terms in the 
expansion close to the critical point vanish. Although the same mechanism is responsible 
for the fact that also even-n representations display non-trivial critical behaviour in the 
d = oo theory, that phenomenon is of course far more difficult to understand from the point 
of view of physics. The even-n Polyakov-line representations simply ought not to be order 
parameters for the deconfinement transition, even in the d = 00 limit, since such sources 
should be screened both above and below the critical point. The resolution of this apparent 
paradox lies in the fact that the critical coupling J c actually vanishes (like 1/d) when 
d—>oo, as follows directly from the mean-field solution (4).f] In terms of the gauge coupling 
g this entails, for fixed N T , g —>■ 00. Although this makes the strong-coupling effective 
Lagrangian analysis more and more accurate, it also pushes the confinement /deconfinement 
phase transition right to the extreme limit g = oo where all sources are "confined" ((Tr n W) = 
for all n at g = 00 in the full gauge theory simply as a consequence of the orthogonality 
property of the group characters). It is for this simple reason that the mean-field solution, 
correctly, predicts critical behaviour for all representations of SU(2). 

The limit d = 00 of finite-temperature gauge theories is thus in many respects highly 
singular. This, together with the Monte Carlo data presented above for the d = 4 SU(2) 
theory, indicates that the usual assumption of d = 00 exponents being valid down to the upper 
critical dimension d u simply fails in this case. Can we understand the singular nature of the 
d = oo limit in an analytical way? As explained above, there are actually no reasons to doubt 
that mean field theory predicts the d = oo behaviour correctly. The only resolution would 
then be that any finite dimensionality d should correspond to radically different behaviour 
close to the critical point, i.e., that 1/rf-corrections discontinuously should alter the critical 
indices. To see whether this is the case, we have considered a slightly improved mean-field 
solution of the same effective Polyakov-line action (1). (This improvement appears to be 
equivalent to what is known as the Bethe-approximation, see, e.g., ref. [p7p. Consider a 



system that consists of 2d +1 sites, a central site (C) and its 2d nearest neighbours (O). 

6 This behaviour is not an artifact of the mean-field solution. It can be checked to hold as well in the 
exact solution of the N — 00 theory p8[. 
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The remaining 2d— I sites of the O-sites are replaced by an external field W, which here is 
the analogue of the conventional mean field. Hence the partition function of this system is 
given by 

Z = [ +1 d$ c 0Ti£n l +1 d$ \Jl-& exp(4J($c + W)$ ) (16) 

J — 1 /-i J — 1 



o 

The integration over the 4> fields leads to 

r+l 



Z oc / d<$>c\ll-<$>l 



hjAJ^c + W)) 12 ' 1 



(17) 



The remaining one-dimensional integration we have performed numerically. In order to fix 
the external field W we require that the magnetization of the fundamental representation is 
equal for the central site (C) and its neighbours (O). We solved this condition numerically 
using the Newton method. Expectation values are evaluated on the central site. 

One can readily check that this improved mean field theory coincides with the conven- 
tional mean field theory in the limit d — > oo. It is, however, expected to be more accurate for 
finite values of d, especially for non-universal quantities. We have solved the self-consistency 
equations for d — 4, 8, 16 and 32. The value found for the critical coupling J c = 0.5352... in 
A-d deviates from the Monte Carlo value by only 2.8%, while standard mean field theory is 
off by 9.2%. But a more striking consequence of the improvement is seen in the behaviour 
of the even-n representations, which are now non-vanishing for all values of the coupling 
J 7^ 0. The numerical values at or below J c , however, decrease rapidly with d. For the 
adjoint representation it is reduced by a factor of approximately 2 at J c when one doubles 
the dimension, while for the n = 4 representation it drops by almost a factor of 4. In this 
fashion the present solution matches the usual mean-field results in the limit d—oo. 

Since such an apparently minor improvement of the mean-field method produces this 
drastic change in behaviour for the even-n representations, it is worthwhile to understand it 
better. In fact, almost any improvement of the d — oo mean-field results are bound to give 
qualitatively the right behaviour for the even-n representations. We will here show that, 
except for the mean-field limit d = oo, the expectation values of Polyakov lines in even-n 
representations will always be positive for all finite d and J. Note that we can rewrite the 
expectation value (Tr n W) as 

(Tr n W) = J dSP(S)M n (S) (18) 

where M n (S) is the conditional expectation value of Tr n W at a site C for fixed sum S of 
the fields <3>o at the neighboring sites O. P(S) is the probability distribution of the sum of 
the fields <E> . M n (S) is given by 



d<S> c exp 


AJSQc + V[<Z> 2 C ] 


Tr 






'AJS^c + V[§ 2 C ] 





MJS) = L - — f r- J — i (19) 

d$> c exp [4 JS<&c + V[<S>1\ 

We have already computed M n (S), since it is of precisely the same form as the mean- field 
solution of the model: we just have to take a = 4JS as argument in eq. (5). For n even 
M n (S) is a even function, since odd (modified) Bessel-functions are odd functions. The only 
real zero of the function is given at S — 0. For all other real arguments M n (S) is strictly 
positive. 
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Next consider the conditional probability distributions p(§no,S), where §no are the 
fields on the neighbouring sites of the O sites. Note that the computation of p(§no,S) 
requires only the integration of the $o fields. One reads from the definition that p($ no, S) > 
for \S\ < 2d for any configuration of the §no with \§no\ < 1- Hence we also have P(S) > 
for \S\ < 2d. 

Taking the properties of M n (S) with n even and P(S) as discussed above (Tr n W) for n 
even has to be strictly larger than zero. Mean-field theory gives zero for J < J c since P(S) 
is replaced by S(S). Any improvement in the mean-field solution that provides a smooth 
function for P(S) will remove the mean-field pathologies that lead to vanishing expectation 
values for the even-n representations at and below J c . Of course, in the limiting case d = oo, 
the exact P(S) is a 5-function at J = J C , and the naive mean-field results are exact. 

With the improved mean field theory we can finally make a much more accurate com- 
parison with our d = 4 Monte Carlo results. In figs. 4 and 5 we have thus plotted (as smooth 
curves) the corresponding predictions for the J-dependent effective exponents defined 
as in eq. (8). Qualitatively the behaviour of the Monte Carlo data is quite well reproduced. 
Since it is clear that the conventional mean field results must be reproduced as d — > oo, 
it is interesting to see what happens with these /^■'■'-exponents as d is increased. In fig. 
6 we show for n = 2 how the lowest-order mean-field behaviour is recovered as d grows. 
One sees that the window close to J c where /3| turns over and starts deviating from the 
d = oo result = \ gets more and more narrow as d is increased. This behaviour is not 
restricted to the adjoint representation. In fig. 7 we show the results for all up to n = 5 
in d — 32 dimensions. The linear spacing of effective exponents is seen throughout, except 
for an extremely narrow interval close to J c . 

Clearly, as d — * oo the window in which the conventional results are reproduced shrinks 
to zero. In fact, one can easily estimate from the improved mean-field solution (16) that this 
window decreases in size as 1/d, eventually disappearing at d = oo. In the more conventional 
language, this is the point at which the amplitudes of the leading terms in the expansion 
for the Polyakov lines vanish. 
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Table 1: Magnetizations at the critical coupling 4J C = 0.5507 for various representations. 
The number in the first bracket gives the statistical error at the given coupling, while the 
number in the second gives the uncertainty due to the error of the critical coupling. 



L 


(T ri W) 


(Tr 2 W) 


(Tr 3 W) 


(Tr 4 W) 


(Tr 5 W) 


4 
6 
8 
12 
16 


0.354(2)(1) 
0.249(2)(2) 
0.194(2)(3) 
0.129(1)(5) 
0.100(1)(7) 


0.1605(14)(6) 

0.1348(7)(6) 

0.1255(4)(7) 

0.1177(2)(6) 

0.1157(1)(7) 


0.0659(4)(2) 
0.0344(3)(2) 
0.0230(2)(3) 
0.0136(1)(4) 
0.0102(1)(6) 


0.0101(6)(1) 

0.00802(30)(6) 

0.00693(14)(6) 

0.00592(8)(6) 

0.00592(4)(6) 


0.0543(3)(0) 

0.02369(14)(1) 

0.01321(7)(1) 

0.00593(3)(1) 

0.00333(2)(1) 



Table 2: Improved mean field theory: J c , and rrij evaluated at J c . 



d 


4J C 


(Tr 2 W) 


(Tr 4 W) 


4 


0.5352319055 


0.07382995665 


0.00254915596 


8 


0.2582840308 


0.03387515014 


0.00055551163 


16 


0.1270085358 


0.01625711255 


0.00013006341 


32 


0.0629953883 


0.00796783664 


0.00003149409 



Figure Captions 

1. ) The fourth-order Binder cumulant for the fundamental representation. The crossing 
determines our best estimate of J c . 

2. ) Same as fig.l, but for the adjoint source. The convergence toward U 2 = 2/3 indicates 
that the expectation value of the Polyakov line in the adjoint representation is non-zero 
throughout. 

3. ) The fourth-order Binder cumulant for the n = 3 source. Although the behaviour is 
radically different from that of the adjoint source, it still does not show the pattern of fig.l 
for the smaller lattices. 

4. ) The effective magnetization exponents (8) for odd-n. The n— 1 and n = 3 representations 
nicely appear to converge toward the Ising value of f3 = 0.5, while the n = 5 representation 
only shows the same trend. The drawn curves refer to the improved mean-field solution 
discussed at the end of the paper. 

5. ) Same as fig. 4, but for the even-n representations. The effective exponent appears to 
converge toward 0, as expected if these magentizations remain non-zero at the critical point. 

6. ) The way the mean-field solution f3 2 = 1.0 is recovered in the limit d^oo. The region in 
{J—J c )/J c where (3 2 eventually turns to zero shrinks as d grows, disappearing in the limit 
d— >oo. 

7. ) The linear spacing of (3^^ , here for <i = 32. 
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